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Abstract: We perform the non-perturbative renormalization of matrix elements of the 
static-light axial current by a computation of its scale dependence in lattice QCD with two 
flavours of massless 0(a) improved Wilson quarks. The regularization independent factor 
that relates any running renormalized matrix element of the axial current in the static 
effective theory to the renormalization group invariant one is evaluated in the Schrodinger 
functional scheme, where in this case we find a significant deviation of the non-perturbative 
running from the perturbative prediction. An important technical ingredient to improve 
the precision of the results consists in the use of modified discretizations of the static 
quark action introduced earlier by our collaboration. As an illustration how to apply the 
renormalization of the static axial current presented here, we connect the bare matrix 
element of the current to the Bg-meson decay constant in the static approximation for one 
value of the lattice spacing, a « 0.08 fm, employing large-volume A'f = 2 data at /3 = 5.3. 
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1. 


Introduction 





In view of the challenging experimental programme of B-factories and the demand of a 
precise quantitative interpretation of its observations within or beyond the Standard Model, 
non-perturbative investigations of the B-meson system and its transition amplitudes in the 
framework of lattice QCD have become a vivid area of research. The impact of lattice QCD 
on this area of flavour physics crucially depends on the precision that lattice computations 
of B-physics matrix elements are able to achieve. Thus, it is very important to try to 
reduce its systematic errors such as the quenched approximation (which currently is already 
being overcome for many phenomenologically relevant quantities, see e.g. [1, 2, 3]) and the 
uncertainties owing to the still unphysically large sea quark masses employed in most 
simulations with dynamical quarks. 

Yet another difficult part of these computations arises from the problem of a sensible 
treatment of b-quarks on the lattice, because lattice spacings small enough to satisfy the 
condition a < 1/mb for a propagating, relativistic b-quark will certainly still continue to be 
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out of reach in the near future. A theoretically clean solution is provided by the heavy quark 
effective theory (HQET). This is a systematic expansion of the QCD amplitudes (between 
hadron states containing a single heavy quark) around the static limit, which describes the 
m\j — ^ DO asymptotics of the effective theory in terms of higher-order corrections multiplied 
by coefficients of 0(l/mb) and powers thereof. For early references and more recent reviews 
consult [4, 5, 6, 7, 8], for instance. Among the attractive features of lattice HQET, the 
theoretically most appealing ones are [9, 10]: (i) Higher-dimensional interaction terms in 
the effective Lagrangian are treated as insertions into static correlation functions, which 
implies that (ii) the continuum limit exists and results are independent of the regularization, 
and (iii) the renormalization of the theory can be performed non-perturbatively, whereby 
also the inclusion of the 0(l/mb)-terms along the basic strategy of [10, 11] has recently 
been implemented in a concrete application [12]. 

However, even the leading (i.e. static) approximation of HQET turns out to be an 
interesting limit, since often it is not expected to be far from results at the physical point 
and, moreover, the static results can yield important information for interpolating in l/mt 
between data at quark masses below the b-quark mass and the static limit. This has been 
demonstrated explicitly for the case of the Bg-meson decay constant, Fq^, in quenched 
QCD in refs. [13, 14, 15], where the value of the decay constant in the static approximation 
was used to constrain the extrapolation of the corresponding heavy-light matrix element at 
finite quark mass values within the charm region. With the present work we want to carry 
out the first step towards a removal of the quenched approximation as one of the main 
systematic errors in the aforementioned determination of Fq^. This step consists in the 
non-perturbative renormalization of the static-light axial vector current in QCD with two 
dynamical quark flavours. In contrast to lattice QCD with relativistic quarks, where the 
renormalization constant of the axial current is only a (lattice spacing dependent) constant, 
the static effective theory gives rise to a scale dependent, multiplicative renormalization 
problem, which can be solved following the strategy of recursive finite-size scaling in an 
intermediate renormalization scheme originally proposed in [16] and already employed for 
the corresponding quenched calculation in [17]. For a review and earlier applications of 
this method, we refer the reader to refs. [8, 18, 19, 20]. 

Let us briefly recall this approach to solve non-perturbative renormalization problems 
in the present context. As for phenomenological applications one is eventually interested in 
matrix elements in a scheme, in which the renormalized amplitudes of the effective theory 
are matched to the QCD ones at finite quark mass, it is usually convenient to first compute 
the scheme independent renormalization group invariant (RGI) matrix element: 



Here, ^*barc = (PS|Aq*^*|0) denotes an unrenormalized matrix element of the static-light 
axial current between a pseudoscalar state and the vacuum, and the renormalization fac- 
tor Zrgi is such that it turns any bare matrix element of ^q*^* into the RGI one. The 
pseudoscalar decay constant at finite quark mass is then related to $rgi through 



$RGI = ^RGl(PS|ArMO) . 



(1.1) 



FpsV^ = Cps {M/A^) X $RGi + O (1/M) , 



(1.2) 
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where M is the RGI mass of the heavy quark and the QCD A-parameter in the 
MS scheme. The function Cps in eq. (1.2) accounts for the fact that in order to extract 
predictions for QCD from results computed in the effective theory, its matrix elements are 
to be linked to the corresponding QCD matrix elements at finite values of the quark mass. 
In this sense Cps translates to the 'matching scheme' [17], which is defined by the condition 
that matrix elements in the (static) effective theory, renormalizcd in this scheme and at 
scale ^ = mb, equal those in QCD up to 1/mb-corrections. Thanks to the three-loop 
calculation of the anomalous dimension of the static axial current in the MS scheme [21], 
the function Cps(M/Aj^) = i*psV"^Ps/^RGi is known perturbatively up to and including 
^^(mb)-corrections. Therefore, the remaining perturbative uncertainty induced in (1.2) by 
the conversion factor Cps is already below the 2% level beyond the charm threshold and 
thus very small. 

Assuming that the running matrix element ^(/u) has been non-perturbatively defined 
in an intermediate renormalization scheme where ji = 1/Lmax represents a low-energy scale, 
the total renormalization factor .^rgi in eq. (1.1) splits into a universal and regularization 
dependent factor according to 



X Zf\gQ,L/a)\^^^^^ X $bare(5o) 

M=l/-^max 



$(//) 

= %Gl(5'o) ^bare(5o) • (1-3) 

The computation of .^rgi(5'o) is the main goal of this work. 

The paper is organized as follows. In section 2 we present our intermediate renormal- 
ization scheme, formulated in terms of the QCD Schrodinger functional. Section 3 contains 
the numerical determination of the (lattice formulation independent) scale dependence of 
the current in this scheme, which is the key prerequisite in order to relate the current 
renormalized at some proper low scale to the RGI current, while section 4 gives our results 
for the (lattice formulation dependent) values of the Z-factor at this low scale. In section 5 
we then explain the use of our results and, as a further illustration, combine them with 
iVf = 2 data for the bare matrix element of the static axial current from ongoing large- 
volume simulations [22] to extract F^^^ for one value of the lattice spacing, a k, 0.08 fm. 
We conclude with a discussion of the results in section 6. Some technical details and tables 
with parts of the numerical results are deferred to appendices. 



2. The renormalization scheme 

We consider QCD with two flavours of mass-degenerate dynamical sea quarks, and the 
heavy quark field is treated in leading order of HQET (static approximation). The renor- 
malization pattern of an arbitrary matrix element $bare of the heavy-light axial vector 
current, 

At^\x) = Ux)l^l^Mx) . (2.1) 

is characterized by the fact that — owing to the absence of the axial Ward identity which 
holds for the corresponding relativistic current — the static-light axial current picks up a 
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scale dependence upon renormalization. Consequently, the scale evolution of the renormal- 
ized matrix element 



= ( PS I {Arnnifi) I ) = PS I V'lToTs^h | ) (2.2) 

between a static-light pseudoscalar state and the vacuum is governed by the renormalization 
group equation 

f^^ = im (2.3) 

d/U 

in formally the same way as it is encountered in conjunction with the running of the 
renormalized quark masses in QCD. 

In the simple form of the renormalization group equation (2.3) we have already implic- 
itly assumed that a mass-independent renormalization scheme is chosen, which is equiva- 
lent to the prescription of imposing renormalization conditions at zero quark mass [23]. 
Moreover, when introducing the lattice spacing a as the regulator of the theory, the 
renormalization factor in question becomes a function of the bare coupling go and a/x, 
^stat _ z^^^(gQ,afi), and only in renormalized quantities this regulator can be removed 
by taking the continuum limit a ^ to finally obtain finite results. In the same way as 
for the renormalized quark masses, also in the static effective theory considered here the 
crucial advantage of mass-independent renormalization schemes is that in all such schemes 
the ratios of renormalized matrix elements constructed as in eq. (2.2) but with a different 
flavour content are scale and scheme independent constants. 

The anomalous dimension associated with the renormalization scale dependence of the 
static-light axial current is encoded in the renormalization group function 7 appearing in 
eq. (2.3), the perturbative expansion of which reads 

7(5) 5' I 70 + 715' +725' + 0(5') I , (2.4) 

with a universal, scheme independent coefficient [24, 25] 

^0 = -4^ (2-5) 

and higher-order ones 71,72, •• • depending on the chosen renormalization scheme. Note, 
however, that generically the 7-function is non-perturbatively defined as long as this is 
the case also for the matrix element $ of the current as well as for the renormalized gauge 
coupling g itself. 

As already emphasized in section 1, we advocate a strategy that regards the renor- 
malization group invariants (RGIs) as the essential physical objects of interest, because 
these are the quantities whose total dependence on the renormalization scale /U vanishes. 
For the present investigation, the relevant RGIs are the renormalization group invariant 
counterpart (1.1) of the matrix element (2.2) of the static-light axial current, 

*ao. = .(.)[2Ml-"-exp{-fd42|-^]}. (2.6) 
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and the QCD A-parameter 



(2.7) 



where the renormahzed coupling g = g{ij,) obeys 

dg 

9- 

with scheme independent one- and two-loop coefficients 



Pig) = fi^ ~ -f\bo + bif + b2g^ + 0{f)} (2. 



Both, $RGi as well as A are defined independent of perturbation theory, and particularly 
the former is not only scale but also scheme independent. 

2.1 Static-light correlation functions in the Schrodinger functional 

A convenient mass-independent renormalization scheme, which has already proven to be a 
theoretically attractive as well as numerically efficient framework to solve renormalization 
problems in QCD similar to that studied here [17, 18, 19, 20, 26], is provided by the QCD 
Schrodinger functional (SF) [27, 28, 29]. It is defined through the partition function of 
QCD on a T X L'^ cylinder in Euclidean space, Z = /^^^3 D[;7, V, V] e-^[^'^''^l, where in the 
lattice regularized form we integrate over SU(3) gauge fields U with the Wilson action and 
two flavours of 0(a) improved Wilson quarks, tp, ip. At times = 0, T Dirichlet boundary 
conditions are imposed on the gluon and quark fields, whereas in the spatial directions of 
length L the fields satisfy periodic (for the quark fields only up to a global phase 9 [30]) 
boundary conditions. Particularly the Dirichlet boundary conditions in time qualify the SF 
as a mass independent renormalization scheme, since they introduce an infrared cutoff to 
the frequency spectrum of qiiarks and gluons and hence allow to perform simulations at zero 
quark mass. The settings T = L, 6 = 0.5 and vanishing boundary gauge fields at xq = 0,T 
then complete the specification of our (intermediate) finite- volume renormalization scheme, 
in which the running renormalization scale /it is now identified with 1/L in a natural way. 

For the non-perturbative renormalization of the static-light axial current in two-flavour 
QCD we closely follow the quenched calculation detailed in [17]. Here, we only recall the 
definition of the basic correlation functions in continuum notation, 

(^o) = -\l d'yd'z Ch(y)75Ci(z)> , (2.10) 

flt\xo) = -\j dVd^z (Mr (^) Ch(y)75Ci(z)> , (2.11) 

/r = / d^ud^vdVd^z <Ci'(u)75Ch'(v) Ch(y)75Ci(z)> , (2.12) 

in terms of the static current (2.1), its 0(a) correction 

5At\x) = (V, + V*) M^) (2.13) 
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Figure 1: Illustration of the correlation functions fx'^^i^n) (left) and (right), defined within 
the Schrodinger functional. The curly and straight lines represent light (i.e. relativistic) and static 
quark propagators, respectively. In the left diagram, the static current Aq*^* is understood to be 
inserted in the bulk of the SF cylinder, at the point where the two quark lines meet. For /|a*(^o)5 
^o*''* is replaced by (JAg**^*. 



and the 'boundary quark and antiquark fields', C, the proper definition of which can be 
inferred e.g. from refs. [31, 32]. These correlators are schematically depicted in figure 1. 

/^^*(xo), for instance, can be shown to be proportional to the matrix element of the (static) 
axial current that is inserted at time distance xq from a pseudoscalar boundary state, while 
the state at the other boundary has vacuum quantum numbers [33]. For the corresponding 
formulae in the lattice regularized theory, in which eqs. (2.10) - (2.12) and any expression 
derived therefrom receive a precise meaning, we again refer to refs. [17, 32]. 

2.2 Normalization condition for the static axial current 

As renormalization condition for the static axial current we impose the condition originally 
formulated in the perturbative context of [32] and later also explored within the non- 
perturbative computation of [17] in the quenched approximation of QCD. Switching to the 
lattice notation from now on, it reads 

Zt\g^,L/a)Xi{gQ,L/a) = Xi{0,L/a) , (2.14) 

where Xi{go, L/a) denotes a suitable 0(a) improved ratio composed from the correlation 
functions (2.10) - (2.12): 

Xi{go, L/a) = --^== . (2.15) 

In this construction, the boundary-to-boundary correlator ff^^^ serves to cancel the un- 
known wave function renormalization factors of the boundary quark fields as well as the 
linearly divergent mass counterterm Sm that one is usually faced with in the static effective 
theory. The definition of Z^'^^ via the condition (2.14) is such that Z^^^ = 1 at tree-level of 
perturbation theory^ (i.e. Zf^^^{0,L/a) = 1); restricting the discussion to the Eichten-Hill 

^The numerical values for the tree-level normalization constant Xi{0, L/a) can be found in table 4 of 
Appendix A in ref. [17]. 
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action for the static quark [5] for the moment, the improvement coefficient cj*^* takes the 
one-loop perturbative value [34, 35] 

cr = -1-00(1) i^5o'- (2.16) 

Moreover, to comply with the demand of setting up a mass independent scheme, 
eq. (2.14) is to be supplemented by the condition of vanishing light quark mass, 

mo = mc , (2-17) 

with bare and critical quark masses as defined in [32]. The critical quark mass fixes the 
hopping parameter, at which the normalization condition for the static axial current at 
given L/a and /3 = Q/g^ is to be evaluated. Its numerical values were already determined 
in ref. [19] from the non-perturbatively 0(a) improved PCAC mass in the light quark 
sector (the latter being defined for 6 = 0.5 through an appropriate combination of light- 
light correlation functions calculated at xq = L/2) and have also been used before in the 
computation of the running of the quark mass in two-flavour QCD [20]. 

Prom the normalization condition (2.14), together with the one-to-one correspondence 
between pairs {L/a,f3), pertaining to a certain fixed value of the renormalized SF gauge 
coupling g^(L), and the box size L as the only physical scale in the system, it is obvious 
that the rcnormalization constant Z^^* runs with the scale jjL = 1/L. Therefore, the change 
of the matrix element $ in eq. (2.2), renormalized in the SF scheme, under finite changes 
of the renormalization scale can now be non-perturbatively computed by means of an 
associated step scaling function o"^^*, which is defined by the change induced by a scale 
factor of two, viz. 

in the continuum limit, it only depends on the renormalized SF coupling u = g^{L). 

The (non- universal) two- loop coefficient of the anomalous dimension (2.4) of the static 
axial current renormalized in the particular SF scheme specified in this subsection is known 
from ref. [32] to be 

71 = if = ( 0.08(2) - 0.0466(13)iVf ) (2.19) 

and will enter the numerical evaluation of the RGI matrix element (2.6) later. 

Before we come to explain the lattice computation of cjJ^*(u) and the subsequent 
steps to relate a bare matrix element of the current to the RGI one, let us comment on 
a difference to the quenched computation of [17]. There it turned out that in case of the 
usual Eichten-Hill action for the static quark the lattice step scaling function Ej'^*(u, a/L) 
(cf. eq. (3.1)) extracted from Monte Carlo simulations acquires large statistical errors in 
the relevant coupling range of ii > 1.5 and that these even grow drastically with L/a. This 
fact originates from the noise-to-signal ratio of the boundary-to-boundary correlator ff^^^ 
in the rcnormalization condition (2.14), which roughly behaves as exp|e*-^)(7Q x (r/a)| 
due to the self-energy of a static quark propagating over a distance T = L. Here, e^^^ 
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is the leading coefficient of the linearly divergent binding energy i^stat of the static-light 
system, .Estat ~ ^ e*-^^ gQ + . . . , and one infers that the precision problem of T,^^^^ becomes 
more severe towards the continuum limit, particularly for the Eichten-Hill action having a 
rather large coefficient, e^^^ = x 19.95 [36]. Thus, for the quenched non-perturbative 
computation, the scheme specified above was finally discarded in favor of a slightly adapted 
scheme [17], in which ff^^^ is replaced by a product of boundary-to-boundary correlators 
involving a light and a static quark-antiquark pair, respectively, whereby especially the 
latter could be calculated with small statistical errors by applying the variance reduction 
method of ref. [37] that consists in estimating the arising one-link integrals by a multi-hit 
procedure. In the case of QCD with dynamical quarks, however, multi-hit can not be 
used, since it does not yield an unbiased estimator any more and — being a stochastical 
procedure rather than analytically defined — it can not be traded for a change in the 
discretization of the action for the static quark. 

Instead of rccoursing to an alternative combination of correlators as done in the 
quenched case [17], we therefore have to pursue a different direction in order to overcome 
the exponential degradation of the signal-to-noise ratio of static-light correlation functions 
computed with the Eichten-Hill lattice action while maintaining the original rcnormal- 
ization condition, eq. (2.14). Fortunately enough, this is indeed possible thanks to the 
alternative discretizations of the static theory devised in rcfs. [13, 38], which lead to a 
substantial gain in numerical precision of B-meson correlation functions in lattice HQET. 

3. The running of the renormaUzed static cixial current 

As emphasized at the end of the previous section, our lattice calculations employ alternative 
discretizations of the static theory that retain the 0(a) improvement properties of the 
Eichten-Hill action [5] but entail a large reduction of the statistical fluctuations of heavy- 
light correlation functions with B-meson quantum numbers [13, 38]. In the following, we 
present results from the static actions denoted as S\^, S^^^ and S^^"^, or for short, 's', 
HYPl and HYP2, respectively. The form and a few properties of these lattice actions are 
briefly summarized in appendix A. 

The light quark sector is represented by non-perturbatively 0(a) improved dynamical 
Wilson quarks, and we refer to ref. [20] for any unexplained details. In particular, the 
dynamical fermion conflgurations, which were generated in the context of that reference for 
a series of given renormalized SF couplings g^{L) at various lattice resolutions, constitute 
the basis for the numerical evaluation of the SF correlation functions (2.10) - (2.12) and 
the renormalization constant Zf^^^ via eq. (2.14). 

Note that Z^^^^ = Zf^^{gQ, L/a) now carries a dependence on the type of static quark 
action used. This holds true also for the step scaling function deduced from it, unless an 
extrapolation to the continuum limit is eventually performed (universality). 
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Figure 2: The lattice step scaling function E^^'(u,a/L) and its continuum limit extrapolations to 
a constant omitting the L/a = 6 data, separately for the three different static discretizations S\^, 
S^^^^ and 5^^^^ (colored lines). The black points (slightly shifted to the left of a = 0) refer to 
the corresponding constrained ^^-minimization of the data from all three actions and give our final 
results for (t^^^{u) in table 1. For the third smallest coupling, the two sets of points at L/a = 6 
refer to one- loop and two-loop values for the boundary improvement coefficient Ct [39, 40]. 



3.1 Continuum extrapolation of the step scaling function 

With eq. (2.18) of section 2 it was already anticipated that the running of a renor- 
malized matrix element of the static axial current in the SF scheme with Nf = 2 massless 
quark flavours can be extracted from the step scaling function a^^^{u), which is defined as 
the continuum limit of the lattice step scaling function S^^*('u, a/L), i.e. 



ar\u) = limEr(n,a/L) , Er(n,a/L) = f Itf f " ' '/^/"^ 
a^o [go, L/a) 



■ (3.1) 
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The condition ttiq = rUc, eq. (2.17), refers to lattice size L/a and amounts to set the hopping 
parameter in the simulations to its critical value, k = Kc- Enforcing g^{L) to take some 
prescribed value u fixes the bare coupling value = 6/P to be used for given L/a. In 
this way, S^^* becomes a function of the renormalized coupling u up to cutoff effects and 
approaches its continuum limit as a/L — > for fixed u. 

We computed Z^'^^{n, a/L) at six values of u, where 
the corresponding box sizes cover a range of the order 
L = 0.01 fm . . . 1 fm (or in 20 GeV ... 0.2 GeV).^ At 
each u, three lattice resolutions — L/a = 6, 8, 12 — were 
simulated, and the numerical results for Z^^^ and Ej^* 
stemming from the three static discretization s, HYPl 
and HYP2 are collected in tables 5 and 6 in appendix B. 
Statistical errors were estimated by a jackknife analysis 
and cross-checked with the method of ref. [41]. 

Given the available data sets belonging to the static 
actions i G {s, HYPl, HYP2} as well as the fact that we 
work in (modulo cj^*, see appendix A) non-pcrturba- 
tively 0(a) improved QCD, there are basically two dif- 
ferent ways to extrapolate the lattice step scaling func- 
tion to the continuum. Either one can perform separate 
fits 

Effiu, a/L) = atf{u) + p,{u) x {a/Lf (3.2) 

or, assuming universality of the continuum limit, a simultaneous extrapolation under the 
constraint of a common fit parameter cr^^* in 

^tf{u, a/L) = at\u) + p',{u) x {a/Lf . (3.3) 

Fit results from two examples of fits of the first kind are deferred to tables 7 and 8 in 
appendix B. These fits employ the ansatz (3.2), including the data at all available lattice 
resolutions in one case, whereas in the other case the coarsest data, L/a = 6, are discarded 
and the fit ansatz is just a constant (p = 0). The nice consistency of the two sets of results 
provides already clear evidence for the very weak overall lattice spacing dependence of 
the step scaling functions Sj^*(u, a/L), particularly beyond L/a = 6, which can also be 
inferred from tables 5 and 6 by direct inspection of the raw data. 

Our final extrapolation to the continuum is based on our experiences gained from the 
quantitative investigations of the running of the QCD coupling and quark masses in the SF 
scheme [19, 20], where the cutoff effects of the corresponding step scaling functions were 
found to be very small as well. The continuum limits were thus obtained according to 
eq. (3.3) by fitting the two values of Sj^* on the finer lattices (L/a = 8, 12) simultaneously 
for i G {s, HYPl, HYP2} to a common constant aj^* (p' = 0), separately for each coupling 
u. We then added linearly the difference between the fit and the L/a = 8 result as a 
systematic error. These constrained constant fits are displayed in figure 2, and the resulting 
^Notc that the exact value of the scale in physical units docs not affect the determination of the renor- 
malization factors in question in this work. As will become clear in the following subsections, it is enough 
to specify the value of the SF coupling for a certain maximal box size in the hadronic regime. 



u 


(Tj'^*(n) 


0.9793 


0.9823(7)(2) 


1.1814 


0.9784(8)(11) 


1.5031 


0.9707(12) (3) 


2.0142 


0.9540(14) (9) 


2.4792 


0.9400(20) (4) 


3.3310 


0.9060(33) (25) 



Table 1: Results for the continuum 

step scaling function a^^^*'{u) from 
constrained fits to a constant, ex- 
cluding the L/a = & data. The first 
error is statistical, while the second 
one is the difference between the fit 
and the L/a = 8 result and is ac- 
counted for as a systematic error. 
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pairs of u and continuum values a^'^^{u) are summarized in table 1. By explicitly trying 
various different extrapolations along the ansatze (3.2) and (3.3) (while not only quadratic 
but also linear in a) we verified with confidence that, within the achievable statistical 
accuracy, our data do not show any significant dependence on the lattice spacing and that 
the continuum limit is well under control. 

3.2 Non-perturbative scale evolution 

For the next steps of the analysis it is most convenient to represent the continuum step 
scaling function a^^^ as a smooth function of u. To this end, we start from eq. (2.6) 
to write down the expression that relates a^^^{u) to the anomalous dimension 7 and the 
/3-function, namely 



/a{u) 



.rM^e.p|-y^ d,Ig.|, (3.4) 

where cr{u) is the step scaling function of the coupling, determined by [19] 

-21n(2)= / dx . (3.5) 

The first non-universal coefficients in the perturbative expansions of 7 and /? are given for 
the SF scheme by the two-loop result 71 as already quoted in (2.19) and the three-loop 
expression (cf. [19]) 



i9)\ 



(4vr) 



^ 0.483 - 0.275A^f + O.OmiN^ - 0.00175A^f^ ) . (3.6) 
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These formulae imply that crj'^* has a perturbative expansion of the form 

af^\u) = 1 + sou + siu^ + S2U^ + ... (3.7) 

with, for instance, the two leading coefficients found to be 

so = ln(2)7o , (3.8) 
si = 2^0 + so ln(2)6o + ln(2)7i . (3.9) 

Guided by eq. (3.7), we therefore represent the non-perturbative results for a^^^^{u) in 
table 1 (with added statistical and systematic errors) by an interpolating fit to an ansatz 
polynomial in u, where so,si are restricted to their perturbative predictions and up to 
three additional free fit parameters are allowed for. All of these fits describe the data 
well, and wc quote the two-parameter fit (the curve of which is shown in figure 3) as the 
final representation for the functional form of a^^^{u). The stability of the polynomial fits 
was further checked by fits with only so (or even no coefficient at all) constrained to its 
perturbative value, which led to compatible results, including a fit to (c7j'^*(u) — l)/n = 
So + siu + ... to reproduce the perturbative prediction for so. 

Moreover, figure 3 demonstrates — in contrast to previous calculations of the non- 
perturbativc scale evolution within the SF scheme of other renormalizcd observables in 
quenched and two-fiavor QCD [17, 18, 19, 20] — that the present case of the = 2 static 
axial current provides another^ example for a significant deviation of the non-perturbative 
data from the perturbative behaviour, which even sets in already at moderate couplings u. 

Now we use (7^^*(u) given by the fit function in order to solve the following joint 
recursion to evolve the coupling and the renormalized matrix element $ from a low-energy 
scale 1/i^max implicitly defined by 

U0=5'(^max) = 4.61 (3.10) 

to the higher energy scales 1/Lk, A; = 0, 1, . . . , 8 (with Lq = Lmax): 

uo = 4.61 , a{uk+i) = Uk =^ Uk = f'{Lk) , = 2"''Lmax , (3.11) 

" k 



For this purpose, also the scale evolution of the coupling is parameterized by an interpolat- 
ing polynomial of the form a{u) = u + aou^ + aiu^ -|- . . . , for which the exact results of the 
corresponding fit (and its covariance matrix) were available from [20, 42]. Since the errors 
on the step scaling functions stem from different simulation runs and are hence uncorre- 
lated, the errors on the fit parameters in the polynomials for a{u),a^^^^{u) and those on 
the recursion coefficients u^, calculated from them can be estimated and passed through 
the recursion straightforwardly by the standard error propagation rules. 



^Comparable deviations between perturbative and non-perturbative running have so far only been ob- 
served for some of the SF schemes studied in the renormalization of four-fermion operators [26]. 
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Figure 4: Numerically computed values of the running matrix element of the static axial current 
in the SF scheme compared to perturbation theory. The dotted and solid lines are obtained from 
eq. (2.6) using the 1/2- and 2/3-loop expressions for the 7- and /3-functions, respectively, as well 
as In(ALmax) = -1.298(58) from ref. [19]. 



3.3 The universal renormalization factor 

Finally, by virtue of eqs. (2.6) and (3.12), we proceed to relate the renormalized matrix 



element ^{fi) - 

$(l/imax) 



'7stat 



-1 



(L)^>barc(5'0) 

Zstat (n—k T \ 

Zstat ( T \ 
A V^tnax ) 



fi = 1/L, at L = Linax to the RGI one as 



(3.13) 



-7o/(2feo) 



exp 



dg 



70 
bog 



with = 2'=/L max and the ratios of Z^— factors, v^, to be taken from the non-perturbative 
solution of the recursion (3.12) discussed in the foregoing subsection. After numerical in- 
tegration of the second factor, ^^Qi/^{1/Lk), with g^ = Uk by employing the perturbative 
expressions for the 7- and /3-functions at two- and three-loop accuracy, respectively, we 
arrive at the series of numbers listed in table 2. For A; > 3 they show a remarkable stability 
in the coupling^ such that we select k = 6 — yielding 



Zstat / r 
A V-'^n 



^A 



0.762(5) at L = 2~^Lr_ 



(3.14) 



^The deviation in the case = is due to the difference between the perturbative and the non- 
perturbative values of ct{u) at large u (see [19]). 
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and the coupling value (uq) entering the second factor in eq. (3.13) to lie safely in the 
perturbative regime — in order to obtain our central result 



£rgi 



0.880(7) or 



^RGI 



1.137(9) at /x=l/Li, 



(3.15) 



$RGl/$(iv 



1 ) 

max/ 



in the SF scheme^ . Hence, through this analysis we have succeeded in connecting the low- 
and high-energy scales L = L^i^^ and L = 2~^Linax non-perturbatively as intended. 

Note that, since the continuum limit has 
been taken, any regularization dependence has 
been removed from the result (3.15) so that 
ideally the error on $(/u)/$rgi of about 0.8% 
should only be included in estimates on ma- 
trix elements in the continuum limit. 

In figure 4 we compare the numerically 
computed running with the corresponding cur- 
ves in perturbation theory. For the argument 
fi/A = l/(LfeA), A; = 0, 1, ... ,8, we plot the 
points $(l/Lfe)/$RGi calculated from (3.13) 
using the universal result (3.15). Here, the 
physical scale A is implicitly determined by 
In(ALinax) = -1.298(58) resulting from the 
recursion (3.11) [19, 20], and the errors of 
the points in figure 4 come from the coeffi- 
cients Vk. As expected from the behaviour 



k 


u 


2/3-loop 


1/2-loop 





4.610 


0.853 


0.851 


1 


3.032 


0.863(3) 


0.862 


2 


2.341 


0.871(5) 


0.870 


3 


1.918 


0.875(6) 


0.874 


4 


1.628 


0.878(6) 


0.877 


5 


1.414 


0.879(7) 


0.878 


6 


1.251 


0.880(7) 


0.879 


7 


1.121 


0.880(7) 


0.879 


8 


1.017 


0.880(7) 


0.880 



of a^^^{u) in figure 3, this comparison be- 



Table 2: Evaluation of eq. (3.13), exploiting 
the perturbative knowledge of the 7- and /3- 
functions in the SF scheme. 



twecn the non-perturbative and perturbative 

running again reveals that good agreement with the perturbative approximation is only 
observed at high scales, while the difference grows up to 5% towards smaller energies of 
the order /x ^ 4A. 



4. Zf^^'^ at the low-energy matching scale 

Connecting a bare matrix element of the static-light axial current to the RGI one according 
to eq. (1.3) amounts to multiply the bare lattice operator with the total renormalization 
factor 

^RGI 



■^RGl(9'o) = 



$(/i) 



X zr(5o,V«)L=^_ , (4.1) 

— l/^/max 



which involves — in addition to the universal ratio $rgi/^(a*)) eq. (3.15) — the value 
of the renormalization factor Zf^^^{go,L/a) at the finite, low-energy renormalization scale 
Lmax implicitly fixed by the condition ^^(Linax) = 4.61 in the intermediate SF scheme. 



^Still, for the perturbative error in this step to be negUgible, it is crucial that 7 is known to two-loop 
precision and /3 to three-loop. 
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L/a 




^stat 


7stat 
^A,HYP1 


ystaX 
^A,HYP2 


„stat 


5.20 


0.13600 


4 


3.65(3) 


0.84621(63) 
0.84846(63) 


0.84550(53) 
0.85760(54) 


0.88820(58) 
0.95061(65) 


1-lp 







6 


4.61(4) 


0.79204(58) 
0.79409(58) 


0.79631(52) 
0.80730(52) 


0.84456(52) 
0.90004(55) 


1-lp 



5.29 


0.13641 


4 


3.394(17) 


0.85323(60) 
0.85541(60) 


0.85218(50) 
0.86385(51) 


0.89202(55) 
0.95227(61) 


1-lp 







6 


4.297(37) 


0.79954(69) 
0.80152(69) 


0.80384(60) 
0.81444(61) 


0.85006(59) 
0.90358(62) 


1-lp 








8 


5.65(9) 


0.75464(71) 
0.75653(71) 


0.76110(65) 
0.77124(65) 


0.80934(65) 
0.86058(67) 


1-lp 




5.40 


0.13669 


4 


3.188(24) 


0.86090(63) 
0.86302(63) 


0.85965(53) 
0.87103(53) 


0.90020(56) 
0.95911(62) 


1-lp 







6 


3.864(34) 


0.81111(68) 

0.81300(68) 


0.81448(60) 
0.82458(60) 


0.85700(59) 

0.90774(63) 


1-lp 








8 


4.747(63) 


0.77310(64) 
0.77491(65) 


0.77816(59) 
0.78785(59) 


0.82316(57) 
0.87185(60) 


1-lp 




Table 3: 


Results for Zf""^ 


with ct set to its 2-loop value, 


both for Ca*"' 


= 1 — loop and cj*^* = 0. 



The values of are from [19]. The hopping parameters k used in the simulations are taken to be 
the critical ones (kc) of [43]. 

Following the steps of the analogous computation in the case of the running quark mass 
in QCD [20], we now derive the second factor in eq. (4.1) for a few values of the lattice 
spacing or the bare coupling, respectively. As pointed out before, this contribution is non- 
universal, and in the form given it will be valid only for our static-light actions, consisting 
of non-perturbatively improved Wilson fcrmions with plaquette gauge action and Cs„ as 
specified in [44] characterizing the light quark sector and the three discretizations s, HYPl 
and HYP2 employed for the static quark flavour. For cj'^* we insert the one-loop values 
recently determined for these static actions in [38] and reproduced in appendix A. It thus 
remains to compute Z^^^^^qq^ L^^la) for the values (3 = 5.2, 5.29, 5.4, which lie well within 
the range of bare couplings commonly used in simulations of two-flavour QCD in physically 
large volumes. The associated simulation parameters and results are summarized in table 3. 
In order to allow for studying the influence of cj*^* on future continuum extrapolations of 
renormalized matrix elements, we also provide estimates for Zf^^^{gQ, L/a) with c^pf^ being 
set to zero. 

While one of the simulations at the largest bare coupling is exactly at the target value 
for (f', the two other series of simulations require a slight interpolation. This has been done 
adopting a fit ansatz motivated by eq. (2.6), 



In (Zr*)=ci+C2ln(52), 



(4.2) 





P 


^stat 


Zrgi 




s 


5.20 


0.7920(6) 

0.7941(6) 


0.6970(6) (56) 
0.6988(6)(56) 


1-lp 






5.29 


0.7873(28) 
0.7892(28) 


0.6928(28)(55) 
0.6945(28)(56) 


1-lp 






5.40 


0.7784(28) 
0.7802(28) 


0.6850(28)(55) 
0.6866(28)(55) 


1-lp 




HYPl 


5.20 


0.7962(5) 
0.8073(5) 


0.7007(5)(56) 
0.7104(5)(57) 


1-lp 






5.29 


0.7922(27) 
0.8026(27) 


0.6971(27)(56) 
0.7063(27) (57) 


1-lp 






5.40 


0.7832(26) 
0.7930(27) 


0.6892(26) (55) 
0.6978(27) (56) 


1-lp 




HYP2 


5.20 


0.8446(5) 
0.9000(5) 


0.7432(5)(59) 
0.7920(5) (63) 


1-lp 






5.29 


0.8390(25) 
0.8919(26) 


0.7383(25) (59) 
0.7849(26) (63) 


1-lp 






5.40 


0.8279(24) 
0.8769(26) 


0.7286(24) (58) 
0.7717(26) (62) 


1-lp 





Table 4: Results for Z^*** and .^rgi for three bare gauge coupling values corresponding to our low- 
energy matching point at g'^ = 4.61 in the SF scheme, distinguishing the three static discretizations 
used. 

to interpolate Z^^^ between those two values of g'^ straddling the target value 4.61, whereby 
the fit takes into account the (independent) errors of both Z^'^^ and g^. We then augmented 
the fit error by the difference between the fit result via eq. (4.2) and the result from a 
simple two-point linear interpolation in g^. The values of the coefficient C2 in the fit (4.2) 
are found to deviate not more than by about 0.03 (with errors on the 7% level) from 
7o/(26o) = -6/29 ^ -0.2069. 

The resulting numbers for Z^^^{go,Linax/a) and finally for the total renormalization 
factor .^RGi(fl'o)) cf. eq. (4.1), are collected in table 4. The first error on Zjiqi stems from 
the error of Zf^^^, whereas the second one embodies the 0.8% uncertainty in the universal 
factor $RGi/^ and, provided that the renormalized matrix element of the static current is 
available at several lattice spacings, should not be added in quadrature to the error on the 
latter before the continuum limit has eventually been taken. For later use, a representation 
of the numerical results of table 4 by interpolating polynomials can be found in table 9 
in appendix B. Comparing the cases c^*^* = 1 — loop and cj*^* = 0, we observe that only for 
the static action HYP2 the change in the renormalization factors is about 6% and thereby 
non-negligible, which however can be attributed to the fact that for this action the one- 
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Figure 5: A^f = 2 conversion factor to the matching scheme, which translates the RGI matrix 
element of Aq*^* to the one at finite mass. A continuous parameterization using the anomalous 
dimension 7 to two- and three-loop accuracy, respectively, is given in eq. (5.4). The dotted curves 
show the corresponding conversion function for A^f = [17, 45] for comparison. 



loop coefficient of cj*^* is by one order of magnitude larger than for the static discretization 
HYPl (cf. eq. (A.5)). 

5. Matrix elements at finite values of the quark mass 

It was already outlined in the introduction that in order to employ results from the static 
effective theory, one has to translate its RGI matrix elements to those in QCD at finite 
values of the heavy quark mass. For the special case of the matrix element of the axial 
current between the vacuum and the heavy-light pseudoscalar, this conversion to the so- 
called ^matching scheme^ amounts to a multiplication with a function Cps(M/Aj^), viz. 

FpsV^ = Cps (M/Aj^) $RGi = Cps (M/Aj^) Z^cii PS | | ) (5.1) 

up to 0(1/M) corrections, where it is theoretically as well as practically advantageous to 
express Cps in terms of a ratio of RGIs as [17, 45] 

n ^A/r/A ^ l7o/(26o) i f^A l-T'^'Hg) Toll 

Cps (M/Aj^) = [ 2bog^ \ exp | d^ ^ -jg^ b^\} ' ^ ^ ^ 

^match (jgnotes the anomalous dimension of the current in the matching scheme. As we will 
see in the next subsection, Cps is very well under control in perturbation theory. 

To exploit eq. (5.1) in order to determine the decay constant Fps, after having non- 
perturbatively solved the renormalization problem of the static axial current, is the main 
purpose of this section. 

5.1 Perturbative evaluation of the conversion function 

The numerical evaluation of the perturbative approximation of the conversion function 
Cps(M/Aj^) has been explained in detail in appendix B of ref. [45]. This analysis is 
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straightforwardly extended to the present case of two-flavour QCD, A^f = 2: Using the 
anomalous dimension in the matching scheme that involves (i) the anomalous dimension 
of the corresponding effective theory operator in the MS scheme up to three-loop order 
[21, 46, 47, 48] and (ii) matching coefficients between the effective theory and physical 
QCD operator up to two loops [5, 36, 49] — together with the four-loop /^-function [50] and 
quark mass anomalous dimension [51, 52] in the MS scheme — , the numerically evaluated 
conversion function Cps(M/Aj^) is shown in figure 5. 

As in the quenched case [17, 45] it is actually more convenient to represent Cps as a 
continuous function in terms of the variable 

X = — 7 : r- (5.3) 

in a form that is motivated by cq. (2.6). This results in 

(3,70/(260) I 1 _ 0.107 a; + 0.093 } 2-loop 7 

, (5.4) 
3.7o/(2bo) I 1 _ 0.118 X - 0.010x2 + 0.043x3 } 3-loop 7 

with bo = (11 — |A^f)/(47r)2 and 70 = — l/(47r2), whereby the prefactor encodes the leading 
asymptotics as x ^ 0. For 3-loop 7, the precision of the parameterization (5.4) is at least 
0.2% for X < 0.6, and one may attribute an error of at most 2% owing to the perturbative 
approximation underlying this determination of the conversion function. 

For future purposes, we also quote the corresponding parameterization for the three- 
flavour theory: 

{a;7o/(26o) { 1 - 0.144X + 0.130x2 } 2-loop 7 

. (5.5) 
a;7o/(26o) { 1 - o.ieix + 0.062x2 + 0.007x3 | 3_iQQp ^ 

5.2 Application: non-perturbative renormalization of F^^* 

For an immediate use of the results obtained, we present a first non-perturbative compu- 
tation of the decay constant Fq^ in the static limit, based on data for the bare matrix 
element of Af^^ in large volume and for the static action S^^^^. To this end we have used 
one of the sets of unqiicnchcd (two degenerate flavours) configurations produced by the 
ALPHA Collaboration for setting the scale for the A^f = 2 A-parameter and RGI quark 
mass [19, 20] by simulations in physically large volumes [22, 53]. The bare parameters are 
given hy 13 = 5.3, k = 0.1355, F = 24^ x 32, and the RGI light quark mass [20] turns out 
to be 171(4) MeV. For the conversion to MeV we have used a = 0.078(1) fm at /? = 5.3 
as given in ref. [54]. The value of the quark mass is at the upper end of the result for Mg 
in [20] (i.e. Mg = 138(31) MeV). Anyway, the dependence of Fb^ on the exact value of 
Mg is expected to be quite mild, as the JLQCD Collaboration, for instance, has reported 
FbJFb^ = l.l3{3){tl^) in ^ef. [55].^ According to that we would have to correct our result 



® A consistent, albeit less precise number was also found in ref. [56] at a finite lattic spacing corresponding 
to (3 = 5.2. 
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on by decreasing it by about 3%. Such a correction is however below the 5% statistical 
error we are able to achieve at this lattice resolution. 

The computational details closely follow what has been done in the quenched case dis- 
cussed in refs. [13, 15]. To suppress excited B-meson state contributions to the correlation 
functions, we introduce wave functions a;(x) at the boundaries of the SF such that the 
correlators Z^^* and /f^* (cf. eqs. (2.10) - (2.12) in section 2) take the form 

fr\xo,co^) = -\ {{Ar%{x)0{u:,)) , fr\u:i,u^) = -\ {O' {co,) 0{u^)) , (5.6) 
with^ 

^M = |lECh(y)75Cc;(y-z)Ci(z), 0'(^) = |!^ci(y)75a;(y-z)caz), (5.7) 

and the improved version (Aj*^*)o of the static-light axial current defined as 

(Af^*)o(x) = A^^^\x) + a4^^SA^^^\x) . (5.8) 
We restrict ourselves to a choice of four spatially periodic wave functions 
'^i(x) = ^ uJi (|x - nL|) , i = l,...,4, 

co.ir) = a-^^/^e-'-/" , a;2(r) = a-^/^e-''/^" , 

a;3(r) = a-^/^e-'^/^a ^ -^^^^ ^ ^-3/2 ^-r/4a ^ (5 9) 

with the coefficients Ni normalizing them such that Q'^ '-^?(-^) ~ -'- bolds. Numerically, 
we have approximated the wave functions using the lowest six Fourier components in each 
spatial (positive and negative) direction. That reduces the computational cost for the 
convolutions required to calculate 

The decay constant is then extracted from the expression for the local RGI matrix 
element of the static axial current, 

$RGi(xo,c.,) = -ZRGi(l + 6ramq)2L3/2-^^e(-o-^/2)^ee(.o,...) , (5.10) 

where for 5j^* the one-loop formula for the static discretization 5'HYP2 £j,qjj^ pgf_ |^3gj enters. 
The effective energy Egs reads 



1 



Eesixo,u}i) = — In 



/^t^t(xo -a,a;i) 



2a [r^''\xo + a,Ui) 



(5.11) 



and is shown in figure 6 for the wave function u>2- This choice is motivated by the fact 
that £'cfj(xo, 0^2) approaches its plateau value earlier than it is the case for the other wave 
functions (which anyway yield consistent values for xq > 1.4 fm). The corresponding bare 
decay constant ^bare; i-e- the quantity in eq. (5.10) with Zrgi set to one and 6a set to zero, 
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Figure 6: Effective energy for the wave function uj2- The plateau region is limited by the dotted 
vertical lines while the horizontal band is the resulting plateau average. 



is displayed in figure 7. Proper linear combinations built from correlators affiliated to the 
other wave functions lead to fully compatible graphs. 

Given the clear plateau in the effective mass, we simply fit it to a constant in the 
region 11 < xq/o < 20 and obtain the result aE'stat = 0.4042(15). Alternatively, we 
perform a two-state fit of aEcf[{xo,u!2) to the function aEgtat + 6i e"^''*''*^" in the range 
^min < xo/a < 20, with tmin > 4. The results for Egtat from the two fits are in complete 
agreement. In addition, in the second case we get numbers for the energy gap aA^^^^, which 
range between 0.5 and 0.7 depending on the choice for tmin- At the same lattice spacing, 
the binding energy aE'stat from the static HYP-actions in the quenched theory turned out 
to be approximately 10% smaller [12, 15]. This is the relative effect one could have guessed 
by using the one-loop expression for the static quark self-energy, considering the shift in f3 
between the quenched and the iVf = 2 theory (i.e. P = 6.1 versus (3 = 5.3, respectively, for 
a ~ 0.08 fm). Thus, the noise-to-signal ratio of the correlation function, which is governed 
by aEsta.t and the mass of the lightest pseudoscalar [58], is comparable in both cases. 

For <I>bare) rather than directly fitting it to a constant within the plateau region, 
we first fit the correlation function ff^^^{xo,u!2) to the two-state ansatz ^2 e~^'''''*^° 4- 
63 e~^'^=''"^^'''''*^^° in order to extract in a second step 'J'bare through the ratio of the coef- 

^In the spirit of ref. [57], in eq. (5.7) we replace one of the spatial sums by a sum over eight separated 
points, which in practice is realized in line with the inversion of the Dirac operator by shifting the source 
at the origin into all octants of the spatial volume L'^. 
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Figure 7: Effective bare decay constant for the same parameters as in figure 6. Note tliat, contrary 
to the effective energy, in this case the horizontal band does not correspond to a plateau average 
but to the value for $barc as obtained from the two-state fits of f^^^{xQ,uj2) (see text). 



ficient 62 and the square root of /f*^*(a;2, ^2), as suggested by formula (5.10). In this way 
rather precise data could be used in the fit, because the main contribution to the statistical 
error of the effective decay constant comes from f^^^^. Moreover, we verified numerically 
that possible excited state contaminations in /^*'^*(a;2, ^2) are strongly suppressed and 
hence can be safely neglected. For varying fit intervals tmin 1^ xo/a < 20 with 4 < tmin < 8, 
this procedure yields stable numbers that are well covered by $bare = 0.1892(62). More de- 
tails will be provided in a forthcoming publication [15] extending ref. [13] which, as stated 
above, we follow quite closely here. 

Combining the result on $bare with eq. (3.15), ^A'^*(fl'o, -^^max/a)|/3=5.3 = 0.8382(25) 
(obtained from the interpolation quoted at the end of appendix B) and am^ = 0.0260(3), 
we get according to eq. (5.10): 

a'/'^RGi = 0.143(5) . (5.12) 

Finally, inserting the proper value of the conversion function Cps(Mb/Aj^) = 1.24(3) (after 
evaluation of the three-loop expression in (5.4) from the previous subsection*^) into eq. (5.1), 
we arrive at F^^^ = 306(14) MeV, where me, = 5368(2) MeV from ref. [59] has been 
employed. The quoted error accounts for the statistical uncertainty and the errors of the 

*In the argument x, eq. (5.3), we take roA^^^^' — 0.62(8) from [19] and tlie recent A'f = result for 
Mb including 0(l/mb)-terms, roMb = 17.12(22) [12]. 
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(universal and discretization dependent parts of the) renormalization factor as well as for 
the errors of the lattice spacing and of Cps- It should be recalled, though, that the light 
quark mass value of the data set analyzed and discussed in this section is still slightly above 
the physical strange quark mass. If the aforementioned decrease by 3% to better meet the 
strange quark mass scale is incorporated, our estimate for the Bg-meson decay constant in 
static approximation at a finite lattice spacing of a 0.08 fm reads 

^stat ^ 297(14) MeV . (5.13) 

This is broadly compatible with the N{ = 2 result F^^^ = 256(45) MeV in the static 
approximation reported in ref. [56] by the UKQCD Collaboration, which however was 
determined at a coarser lattice spacing of a ^ 0.1 fm (/3 = 5.2) and which rests upon the 
Eichten-Hill static quark action and only perturbative renormalization. 

We can get an idea about discretization effects from the quenched (iVf = 0) compu- 
tation [15] of the very same quantity, again using 0(a) improved Wilson fermions and 
plaquette gauge action but the HYPl action for the static quark. There, Fg*^* in the con- 
tinuum turned out to be 233(18) MeV, and in that case the value of F^^^ becomes about 
10% larger at a ~ 0.08 fm than the continuum limit [15]. Therefore, our result indicates an 
increase in the two-flavour theory either for the value of F^^ , as already observed in ref. [3] , 
or for the 0(l/mb)-corrections to the decay constant. 

6. Conclusions 

We have presented a solution of the scale dependent renormalization problem of the static 
axial current in two-flavour QCD by means of a fully non-perturbative computation of 
the renormalization group running of arbitrary renormalized matrix elements $ of Aq'^^ in 
the Schrodinger functional scheme.^ In particular, the renormalization factor $rgi/^(/^) 
relating the matrix element at a specified low-energy scale /i with the associated RG invari- 
ant in the continuum limit, eq. (3.15), is obtained with a good numerical precision, which 
is comparable with the corresponding study in the quenched approximation [17]. This is 
an important prerequisite for a controlled determination of Fb in the static limit of the 
two-flavour theory. 

The use of the alternative discretization schemes of refs. [13, 38] for the static quark 
instead of the traditional Eichten-Hill action has not only led to a substantial reduction of 
the statistical errors of the static-light correlation functions involved in our computation, 
but also entails a convincing universality test of the continuum limit of the current's lattice 
step scaling function (cf. figure 2). Similar to the case of the renormalized quark mass 
in ATf = 2 [20], we find this function to be nearly independent of the lattice spacing for 
a/L < 1/6 and hence conclude that 0(a) improvement is very successful also at large values 
of the SF coupling. Contrary to the running of the quark mass, however, the scale evolution 

® As a consequence of the heavy quark spin symmetry that is exact on the lattice and owing to the chiral 

symmetry of the continuum theory, our computation even yields the scale dependence of all static-light 
bilinears ^'i^i'h, up to a scale independent, relative renormalization [17]. 
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of the non-perturbative step scaling function of the renormahzcd static-hght axial current, 
figure 3, exhibits a pronounced deviation from perturbation theory already at moderate 
couplings. 

As a first application, we have combined our non-perturbative renormalization factors 
with a numerical result for the bare matrix element of ^q*^* extracted from a large- volume 
Nf = 2 simulation, in order to estimate the Bg-meson decay constant in the static limit for 
one value of the lattice spacing, a ~ 0.08 fm. In QCD with dynamical quarks, the static 
approximation may even be the starting point of the most viable approach to determine 
Fb, especially if such a computation is supplemented with data on the heavy-light decay 
constant from the charm quark mass region (as already demonstrated in the quenched 
case [13, 15]) or with 1/M-corrections. For a controlled inclusion of the latter, also the 
matching between the effective theory and QCD inherent in the conversion function Cps 
(cf. eq. (5.1)) should then be performed non-perturbatively [8]. A general framework for 
this is provided by the non-perturbative formulation of HQET exposed in ref. [10]. Other 
approaches, such as the method of heavy quark mass extrapolations of finite- volume effects 
in QCD [60] and its conjunction with HQET [61], have been proven to be feasible in 
quenched QCD and yield consistent results there, but it may turn out to be difficult to 
extended them to the dynamical case. 

With -Fg*^* = 297(14) McV at a 0.08 fm we find a value for the static decay constant 
that is significantly larger than the A^f = estimate F^^ = 194(6) MeV quoted as the 
best quenched result in the recent review by Onogi [3]. This signals a visible effect of the 
dynamical fermions in the two-flavour theory — qualitatively in line with the outcome 
of other unquenched calculations [3] — which in the end could reflect in an increase of 
i^Bs and/or its 1/M-corrcctions to the static limit. Of course, these issues deserve further 
investigations and can only be settled if the two remaining sources of systematic errors 
are overcome, namely the yet slightly too large value of the sea quark mass and the lack 
of a continuum limit for F^^^. Therefore, the reduction of the sea quark mass as well as 
a determination of the bare static-light decay constant for a smaller lattice resolution are 
part of an ongoing project [22]. Rather in the long term, also an active strange quark 
(Af = 2-1-1) will have to be included. 
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A. Lattice actions for the static quark sector 

We use three discretizations of the action for static quarks, introduced in refs. [13, 38] 
to yield an exponential improvement of the signal-to-noise ratio in static-light correlation 
functions compared to the standard Eichten-Hill action [5] , 

5h[f/,?h,V'h] = a4 5]?h(^)^oV'h(a;) , (A.l) 

X 

Dq being the time component of the backward lattice derivative acting on the heavy quark 
field il)\i{x). These new static quark actions rely on changes of the form C/(x, 0) W{x,Q) 
of the parallel transporters U{x, fj,) in the covariant derivative, 

Doil^h{x) = ^ iph{x) - W\x - aO, 0)iph{x - aO) , (A.2) 

where now W(x,0) is a function of the gauge fields in the immediate neighborhood of x 
and x + aO. In the numerical simulations and the data analysis underlying the present work 
we have considered, among the possible choices^*^ for W{x,0), the regularized actions 



^2 

W{x,0) = V{x,0) 



^ + (^lTrVHx,0)V{x,0)y^' 



'h ) 



t^(x,0) =Vhyp(x,0) => ^HYP, (A.4) 



(A.3) 



where V(x,0) is the average of the six staples around the link U{x,0) and VHYp(a;,0) the 
HYP-link, the latter being a function of the gauge links located within a hypercube [62]. 
The 'HYP-smearing' involved in the construction of the HYP-link requires to further specify 
a triplet of coefficients, the two choices {ai, 02, a^) = (0.75,0.6,0.3) and {ai, 02, as) = 
(1.0, 1.0, 0.5) of which were motivated in [62] and [38], respectively, and define the associated 
static actions 5^^^-'^ and S"^^^^. The discretization is inspired by the SU(3) one-link 
integral, for which W{x, 0) in the form of eq. (A.3) is an approximation. For more details 
the reader may consult ref. [38]. In the text, we will also frequently distinguish these three 
static actions by just referring to them as 's', HYPl and HYP2. 

While the largest improvement in statistical precision is actually achieved for the action 
HYP2, it is observed [13, 38] that generally with all the proposed discretizations at least 
an order of magnitude in the signal-to-noise ratios of B-meson correlation functions in the 
static approximation can be gained at time separations around xq 1.5 fm w.r.t. the action 
(A.l) and that, in addition, the statistical errors grow only slowly as xo is increased. Even 
more importantly, quite the same small scaling violations in the 0(a) improved theory are 
encountered with the new discretizations [38]. 



^"The sensible choice of the parallel transporters is guided by the demand of preserving on the lattice 
those symmetries of the static theory, which guarantee that the universality class and the 0(o) improvement 
are unchanged w.r.t. the Eichten-Hill action [13, 38]. 
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In that reference, also the (regularization dependent) improvement coefficient cj*^* 
multiplying the 0(a) correction (2.13) to the static-light axial current has been numer- 
ically determined in one-loop order of perturbation theory, and we here reproduce the 
corresponding expansions for the three static quark discretizations at our disposal: 



cr(5o) 



stat,(l) 



X 5q with 



stat,(l) 



0.0072(4) for 3% 
0.039(4) for 5^ypi 
0.220(14) for^HYPS 



(A.5) 



B. Simulation results for Z^'^^ 

In tables 5 and 6 we collect the bare parameters and results of our simulations to compute 
^|tat rpj^g pairs {L/a,l3) and the associated values of the critical hopping parameter, 
K = Kc, were already known from the non-perturbative computation of the running of the 
Schrodinger functional coupling itself [19]. 

In order to extract the lattice step scaling function Sj*^* according to eq. (3.1), simu- 
lations on lattices with linear extensions L/a and 2L/a are required. At the three lowest 
couplings g'^{L), the runs have been performed using the one- loop value of the boundary 
0(a) improvement coefficient Ct in the gauge sector [39], 



cl-'^{go) = l-Omig'o , 

except for /3 = 7.542 at L/a = 6 and /3 = 7.7206 at L/a 
as for the larger couplings the two-loop value of Ct [40], 



(B.l) 

8. For these parameters as well 



c^'P(5o) = 1-0.05150 - 0.0305^ 



(B.2) 



has been employed. At the third lowest coupling, u ^ 1.5, we checked at L/a = 6 that 
there is no significant difference in E using the one- or two-loop value for ct . This is also 
illustrated by the additional data points included in figure 2. For decreasing a/L, the effect 
of the accuracy in ct on the results is expected to become even smaller. The contribution 
to the error of S^*^* induced by the uncertainty in the coupling u (which can be estimated 
with the aid of the one-loop result ln(2)7o for the derivative of cr^^* with respect to u) is 
negligible compared to the statistical error of S . 

Tables 7 and 8 list the results of the alternative continuum hmit extrapolations of the 
lattice step scaling function mentioned within the discussion in Section 3.1, which were 
performed as separate fits of the three data sets corresponding to the static actions s, 
HYPl and HYP2. 

Finally, we summarize in table 9 the coefficients Zi, fi of the polynomial interpolations 
as functions of 5.2 < /3 = G/g^ < 5.4, 



^stat 



{go, L^s^/a) = J2zi{P- 5.2y , Z^Giigo) = (/? - 5.2)^ 



(B.3) 



i=Q 



1=0 



^^The SF-specific boundary 0(o) improvement coefficient that involves the quark fields, Ct, was set to its 
one- loop value [63] throughout. 
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of the numerical results on the renormalization factors Z^^*(5(o, Lmax/a) and ^rgiIs'o) 
tabulated in section 4. The statistical uncertainty to be taken into account when using 
these formulae varies between 0.1% (/3 = 5.2) and about 0.3% (/3 = 5.4). Only in the case 
of Z^Gi, the additional 0.8% error of its regularization independent part (3.15) needs to be 
included. 



OS 





B 

y 




IJ 1 Li 








0.9793 


9.50000 
9.73410 
10 05755 


0.131532 
0.131305 
131069 


6 
8 
1 9 


0.9396(5) 
0.9321(6) 
Q953f'4'l 


0.9190(7) 
0.9146(9) 
QOSSf?"! 


0.9782(8) 
0.9813(11) 
Q81 fifQ"! 


1.1814 


8.50000 
8.72230 
8 QQ3fifi 


0.132509 
0.132291 
131Q75 


6 
8 
1 9 


0.9277(4) 
0.9195(7) 
Ql 1 ^(A^ 


0.9030(9) 
0.8997(10) 

sQfi'^n 1 


0.9733(10) 
0.9785(13) 
Q7fi8fl 91 


1.5031 


7.50000 
8.02599 


0.133815 
0.133063 


6 
12 


0.9095(6) 
0.8922(8) 


0.8783(14) 
0.8649(14) 


0.9658(17) 
0.9694(19) 


1.5078 


7.54200 
7.72060 


0.133705 
0.133497 


6 
8 


0.9110(6) 
0.9012(11) 


0.8780(9) 
0.8736(13) 


0.9638(12) 
0.9694(18) 


2.0142 


6.60850 
6.82170 

7.09300 


0.135260 
0.134891 

0.134432 


6 
8 
12 


0.8843(9) 
0.8752(15) 

0.8636(11) 


0.8422(11) 
0.8343(16) 
0.8219(17) 


0.9525(15) 
0.9533(24) 
0.9518(23) 


2.4792 


6.13300 
6.32290 
6.63164 


0.136110 
0.135767 
0.135227 


6 
8 
12 


0.8623(11) 
0.8506(16) 
0.8417(9) 


0.8021(18) 
0.7984(25) 
0.7894(24) 


0.9302(24) 
0.9386(34) 
0.9379(32) 


3.3340 


5.62150 
5.80970 
6.11816 


0.136665 
0.136608 
0.136139 


6 
8 
12 


0.8310(13) 
0.8152(14) 
0.8061(14) 


0.7450(22) 
0.7367(36) 
0.7333(35) 


0.8965(29) 
0.9037(47) 
0.9097(46) 



Table 5: Results for the step scaling function S^'^* with discretization S\. 
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Table 6: Results for the step scaling function SJ^*** with discretizations HYPl and HYP2. 
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Table 7: Results of the continuum limit extrapolation of the lattice step scaling function 
Y?^^^{a/ L,u) to C7^'''(u), fitting the data from our three discretizations separately at all available 
lattice resolutions as a function linear in (a/L)^. 
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Table 8: As in table 7 but upon omitting the L/a — 6 data and fitting to a constant. 
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Table 9: Coefficients of the interpolating polynomials of the renormalization factors in eq. (B.3). 
Uncertainties are discussed in the text. 
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